Phase-shifts in stochastic resonance in a Chua circuit 
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We present an experimental study of stochastic resonance in an electronic Chua circuit operating 
in the chaotic regime. We study in detail the switch-phase distribution and the phase-shift between 
sinusoidal forcing for two responses of the circuit: one depending on both inter-well and intra-well 
dynamics and the other depending only on inter-well dynamics. We describe the two relevant de- 
synchronizatrion mechanisms for high and low frequencies of the forcing and present a method to 
detect the optimal noise intensity from switch phases which coincides with the one derived from the 
observation of the signal-to-noise ratio or residence times. 

PACS numbers: 05.40.-a,02.50.Ey 



I. INTRODUCTION 

Stochastic resonance, by which a weak signal acting 
upon a nonlinear system can be amplified by the addition 
of noise of optimal intensity, has attracted considerable 
attention in recent years. While most of the research has 
focused on bistable and excitable systems [l], Q , studies 
have also been performed in chaotic systems with two at- 
tractors, in what can be considered as a generalized form 
of a bistable system 0, 0] ■ Several measures have been 
introduced to quantify the strength of the response to the 
forcing: the response amplitude originally used by Benzi 
et al. HI , the signal-to- noise ratio [if, H[ , as wen as other 
measures determined from the residence and switching 
times probability distributions @, 0, H[ or the work done 
by the external force @. Although all these different 
quantifiers display a maximum as a function of the noise 
intensity, the values of the optimal noise intensity ob- 
tained in each case do not necessarily coincide. 

In both theoretical and experimental studies the phase- 
shift between the force and the response has been much 
less considered, although it is interesting from the view- 
point of relating stochastic resonance to standard res- 
onance phenomena. Moreover, the magnitude of the 
phase-shift and its variation with noise intensity has 
been the subject of some controversy 
phase-shift has been analytically calculated in bistable 
system within the linear response theory taking into ac- 
count both intra-well and inter- well dynamics [l(| EH- 
It displays a bell-shape dependence on the noise inten- 
sity with a extremum at a smaller value of the noise in- 
tensity than the optimal noise intensity obtained from 
the signal-to-noise ratio. This result was confirmed by 
measurements on an analog electronic circuit [l4[ and 
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the same dependence has also been obtained numerically 
[^.[l5|. However, results in a one-dimensional Ising model 
show a monotonous change of the phase-shift with the 
noise intensity fl6ll. This controversy [111, [12] was con- 
sidered in Refs. jl6l. \vK . [H| based on numerical studies 
of a bistable system and an assembly of superparam- 
agnetic particles. One obtains the monotonous change 
of the phase-shift with the noise intensity if only inter- 
well dynamics is taken into account. The bell-shape de- 
pendence thus reflects the competition between inter-well 
and intra-well dynamics. It was then concluded that the 
extrema in the dependence of the phase-shift and signal- 
to-noise ratio have a different origin. The phase-shift has 
also been studied in terms of switch-phase distributions 
in Ref. 19]. The authors presented a numerical investiga- 
tion in the symmetric Schmitt trigger and proposed the 
de-synchronization mechanism responsible for the disap- 
pearance of stochastic resonance. 

The aim of this paper is to present experimental results 
of our thorough studies of the phase-shift in stochastic 
resonance in the Chua electronic circuit [20[ | operating in 
a chaotic regime where two single-scroll attractors coex- 
ist. In this case the circuit can be thought of as bistable 
and the observed phenomena correspond to conventional 
stochastic resonance within two stable wells [1J]. We have 
determined the phase-shift between a sinusoidal forcing 
E(t) = E sin(ujt) and two responses of the circuit: one 
depending on both inter-well and intra-well dynamics and 
the other depending only on inter-well dynamics. This 
clarifies the dependence of the phase-shift on the noise 
intensity and parameters of external forcing. Consider- 
ing switching-phase distributions one can distinguish two 
de-synchronization mechanisms for low and high frequen- 
cies of external forcing. This leads us to propose different 
methods to determine the optimal noise intensity. We 
show that the same optimal noise intensity can be ob- 
tained from the signal-to-noise ratio and the switching, 
residence time or switching-phase distributions. Stochas- 
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FIG. 1: Diagram of electronic Chua's circuit. E(t) is the forc- 
ing sinusoidal signal and £(t) is a Gaussian white noise source. 
R = 16750, L = 18mH , d = lOnF and C 2 = lOOnF . The 
nonlinear Chua's diode Nr has the current- voltage character- 
istic: / = /(V) = mil/ + l(77To - mi){\V + V \ - \V - Vo\), 
with parameters mo = — 0.758mA/V, mi = — 0.409mA/V 
and Vo = 1.08V. Vi is the voltage on the capacitor Ci . 



tic resonance is a widespread phenomenon and we hope 
that our methods will provide an easy and appropriate 
way to determine numerically or experimentally the op- 
timal noise intensity. 



II. EXPERIMENTAL SETUP AND RESULTS 

The Chua circuit and parameters of its components 
used in our experiments are given in Fig[T] Stochastic 
resonance in this circuit was observed and described in 
Ref.[l|- In the absence of the forcing, E(t), and the noise 
source, <f(t), the dynamics displays two single-scroll at- 
tractors and trajectories flow to one or another depending 
on initial conditions. Jumps between the attractors can 
be observed if the amplitude Eq of the forcing is above a 
threshold value that depends on the frequency to of the 
forcing [4]. In all experiments we first set the amplitude 
of forcing just below that threshold value and then add 
the noise signal to the forcing so to induce jumps of the 
system between single-scroll attractors. During the ex- 
periments we recorded both E(t) and the voltage V\(t) 
on the capacitor C\. This voltage is selected because it 
clearly shows the position of the dynamical trajectory 
on either attractor and it represents the response of the 
circuit depending on both inter-well and intra-well dy- 
namics. The step function representation, S(t), of this 
voltage has been used in many studies [l[ as it represents 
the response of the circuit depending only on inter-well 
dynamics. 

We have computed phase-shifts both for V\{t) and 
S(t). An example of temporal evolution of E(t), V\{t) 
and S(t) is shown in FigO At forcing frequencies be- 
low 100Hz one can determine the phase-shift between 
the forcing E{t) and the response, e.g. V\if), as the value 
$Vi °f the phase <f> such that the cross-correlation func- 
tion (E(t)Vi(t + (fr/cj)) takes the maximum value (where 



FIG. 2: Temporal evolution of the forcing E(t) = Eq sm(ut) 
and the voltage Vi(t) with its step function representation 
S(t). The amplitude of the forcing signal is Eq = 535 mV 
and the frequency to — 7.5 Hz. 



(• • •} represents a time average). In Fig. Owe present the 
dependence of the phase-shift $vi on the noise intensity 
(measured as the standard deviation of its probabiity dis- 
tribution) for different amplitudes Eq and frequencies u>. 
This is the bell-shape dependence observed and explained 
by Dykman et al. [Hj]. For small noise levels, inter-well 
hoppings can be neglected and the circuit dynamics de- 
pends almost totally on the interaction between E(t) and 
the intra-well motion. The characteristic frequency of 
the dynamics of our circuit at the single-scroll attractor 
is around 2740Hz, so the phase-shift remains small. The 
abrupt decrease in the phase is associated with the onset 
of inter-well jumps. 

One can notice that the minimum of the phase-shift de- 
pends both on Eq and u. This is illustrated by the data 
presented in Figj3] in the case u> = 1Hz. The threshold 
amplitude value in this case is about 555mV. For lower 
forcing amplitudes higher noise intensities are necessary 
to induce jumps between attractors and the minimum 
of the bell-shape curve moves to higher noise intensities. 
At low u> the phase-shift $5 between E(t) and the step 
function representation of the voltage S(t) can also be 
determined from the maximum of the corresponding cor- 
relation function. An example of the dependence of this 
phase-shift on the noise intensity is shown in FigU in 
the case u> = 10Hz. This dependence is monotonic and 
similar to what has been obtained in a one-dimensional 
Ising model [l6j . By superimposing in this figure also the 
phase-shift $Vi one can notice that as soon as inter-well 
jumps are activated they immediately become dominant 
and cause the decrease of the phase-shift between E(t) 
and V\(t). The inter- well hoppings also make the phase- 
shift tend to zero for higher noise intensities. 

Phase-shifts can also be determined from the switch- 
phase distribution. This method has been used in the 
symmetric Schmitt trigger [l9| and it is specially impor- 
tant for theoretical and experimental studies of stochastic 
resonance in neurons [2l|. In our case we looked at the 
distribution f + {<j>+) of the oriented switch phases </>+ de- 
fined as the phases of the sinusoidal forcing signal (mod- 
ulo 27r) corresponding to switches of the step function 
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FIG. 3: The dependence on the noise intensity of the phase- 
shift $vi between the forcing E(t) and the output voltage 
Vi(t) for different amplitudes and frequencies. 



FIG. 5: Switch-phase distribution for different noise intensi- 
ties (see the legend). The amplitude of the forcing signal is 
535mV and the frequency 7.5 Hz. 




standard deviatidn of noise [mV] 

FIG. 4: Phase-shifts <J?Vi , 5>s (defined from the correlation 
of the forcing E(t) and the responses Vi(t) and S(t), re- 
spectively) and the average switch-phase <&+ plotted as a 
function of noise intensity. The amplitude of the forcing is 
Eq = 505mV and the frequency u = 10Hz. 

S(t) from negative to positive value. The distribution 
/_(<^_) of the switch phases 0_ in the opposite direction 
can be obtained by a translation by n of the distribution 
of the switch phases /+(</>+). The average of the switch 
phases <£> + is defined by the following equation: 

1 N 

pexp (»*+) = — exp (»$.) (1) 
fc=i 

where <K_, with k — l,...,N, are the different phase 
switches observed during the time evolution and the am- 
plitude p of the complex quantity serves as an order pa- 
rameter characterizing the degree of phase synchroniza- 
tion in neurons [2l[ . The dependence of <&+ on the noise 
intensity coincides with that of the phase-shift $5 be- 
tween E{t) and S(i), as shown in FigfJJ Both measures 
can thus be equivalently used to determine the phase- 



shift of a system whose response depends only on inter- 
well dynamics. In Figj5]we present distributions of switch 
phases /+(</>+) obtained at forcing frequency u) = 7.5Hz 
and different noise intensities. At low noise intensities 
there is only one peak in the distribution. As the noise 
intensity increases, this peak shifts towards zero and flat- 
tens. For a standard deviation of noise around 300mV 
a second peak appears at distance n from the domi- 
nating peak. The onset of the second peak signals the 
onset of de-synchronization in stochastic resonance. As 
the noise intensity increases the second peak grows until 
both peaks become comparable. This de-synchronization 
mechanism was observed and described in the numerical 
studies of the symmetric Schmitt trigger [n| and signals 
the degradation of the stochastic resonance with the noise 
intensity. 

These observations allow us to propose a quantity 
which can be used to determine the optimal noise level 
from switch-phase distributions. Let us denote by P<f>, 
the probability that the switch-phase <p + belongs to the 
interval — tt/2,4)® + ir/2], where is the posi- 
tion of the dominating peak in the switch-phase distri- 
bution. The dependence of this probability on the noise 
intensity for different forcing frequencies is presented in 
FigjS] For frequencies higher than 7.5Hz all curves in 
this figure have the inflection point at a noise inten- 
sity of 300mV. In a previous paper we considered 
the residence time distributions and defined an appro- 
priate quantity which is suitable to detect the optimal 
noise level from these distributions. This quantity is the 
probability Pt t that the residence time is in one of inter- 
vals [(« - 3/4)^, (i - 1/4)^], where i = 1,2,.... The 
dependence of this probability on the noise intensity for 
the same forcing frequencies as in Figgis shown in FigjT] 
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FIG. 6: The probability P$ + that the switch-phase $ + is in 
the interval - tt/2, $J + vr/2] as a function of noise inten- 
sity, denotes the position of the dominating peak in the 
switch-phase distribution. The legend gives the amplitudes 
and frequencies of the forcing periodic signal. 
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FIG. 7: The probability Pr r that the residence time is in one 
of intervals [(i - 3/4) 2s, (i - 1/4) 2^] , where i = 1, 2, . . . , as 
a function of noise intensity. The legend gives the amplitudes 
and frequencies of the forcing periodic signal 



All curves at this figure also have the inflection point for 
the same noise level as curves in Fig|5] We observed [4j 
that this noise level is the same as the optimal noise level 
obtained experimentally from the dependence of signal- 
to-noise ratio on noise intensity. We thus propose to use 
alternatively the dependence of the signal-to-noise ratio 
on the noise intensity or the dependence of probabilities 
P$ + or Pt t on the noise intensity. In the first case the 
optimal noise level corresponds to the maximum whereas 
in the other cases it corresponds to the inflection point. 

The methods described above to determine the phase- 
shift and de-synchronization mechanism change for forc- 
ing signals with frequencies above 100Hz. In this case, 



FIG. 8: Switch-phase distribution /+(</>+) for different noise 
intensities. The amplitude of the forcing signal is 14mV and 
the frequency 540Hz. The standard deviations of noise are 
given in the legend. 

the correlation function between the response of the sys- 
tem (either V\(t) or S(t)) and the forcing E(t) vanishes 
and the synchronization between the forcing and the re- 
sponse can only be determined from the switch-phase 
distributions, as the ones presented in Fig[5] for differ- 
ent noise intensities and u> = 540Hz. One can note that 
there is only one peak in the distribution for all noise 
intensities. Increasing the noise intensity the position of 
the peak does not change but the distribution flattens. 
Finally for very large noise intensities the distributions 
of switch phases </>+ and </>_ start to overlap. This is 
the second mechanism of de-synchronization that ensues 
stochastic resonance on increasing further the noise level. 
The position of the peak in the switch-phase distribution 
changes with the forcing frequency. We present in Fig[9] 
the dependence on the forcing frequency of the average 
switch-phase $ + , as obtained from Eq. |T]). The values 
<E>+ = — 7r/8, — 7r/4, — 7t/2 correspond, respectively, to fre- 
quencies u> — 590Hz, 1080Hz, 2160Hz. This dependence 
does not depend on the value of the sub-threshold ampli- 
tude of the forcing signal. It is evident from Figs. Elands 
that also for high frequencies the optimal noise intensity 
is determined by the inflection point of the dependence 
of the probability P r j >+ or Pt t on the noise intensity. As 
we have shown in Ref. [|| this corresponds to the opti- 
mal noise intensity as determined by the maximum of 
the signal-to-noise ratio. 

III. CONCLUSIONS 

In this paper we have addressed the study of the the 
phase-shift in stochastic resonance. We have analyzed 
two mechanism of de-synchronization observed at high 
and low forcing frequencies. Our results point out the 
attention to switch-phases which are able quantify the 
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synchronization between the forcing and response in the 
whole frequency range of the forcing signal. This paper 
also suggests methods to detect the optimal noise inten- 
sity by observing the signal-to-noise ratio, the residence 
times, the switching times or switch-phases. The optimal 
noise intensity determined from any of these observations 
is the same and they can be used alternatively. We have 
derived our main results by studying experimentally a 
Chua circuit operating in a chaotic regime, but we be- 
lieve that our main conclusions are rather general and 
could be useful in many investigations and applications 
of the phenomenon of stochastic resonance. 



FIG. 9: The dependence of the average switch-phase $+ on 
the frequency of the forcing sinusoidal signal with an ampli- 
tude close to the threshold amplitude. 
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